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Abstract 

Recent measurements at RHIC suggest that a nearly perfect fluid of quarks and gluons is pro- 
duced in collisions. Moreover the passage of supersonic partons through this medium seems to 
produce waves. These waves might pile up and form Mach cones, which would manifest themselves 
in the so called away-side jets, forming a broad structure in the angular distribution of the particles 
recoiling against a trigger jet of moderate energy. In most of the theoretical descriptions of these 
phenomena, the hydrodynamic equations are linearized for simplicity. We propose an alternative 
explanation for the observed broadening of the away-side peak. It is based on hydrodynamics but 
it is a consequence of the non-linearities of the equations, which instead of simple waves may lead 
to localized waves or even solitons. 

We investigate in detail the consequences of including the non-linear terms. We use a simple 
equation of state for the QGP and expand the hydrodynamic equations around equilibrium con- 
figurations. The resulting differential equations describe the propagation of perturbations in the 
energy density. We solve them numerically and find that localized perturbations can propagate for 
long distances in the plasma. Under certain conditions our solutions mimick the propagation of 
Korteweg - de Vries solitons. 



I. INTRODUCTION 



The heavy ion colhsions performed at BNL's Relativistic heavy Ion Colhder (RHIC) 
create a hot and dense medium, which behaves as a perfect fluid. During the first years 
of the RHIC program, hydrodynamics was applied to describe the space-time evolution 
of the bulk of the fluid. In the last years hydrodynamics became relevant to study also 
the perturbations on the fluid, such as, for example, the waves generated by the passage 
of a supersonic parton. This field was opened by the observation of a broad structure in 
azimuthal di-hadron correlations P, 9]. This broad structure is called the "away-side jet" 
and recoils against the "near-side jet" (or "trigger jet"). In the framework of hydrodynamics, 



this observation could be explained by the conical shock waves gen erated 



deposition in the hydrodynamical medium 
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jy large energy 
1^. Although 



quite elegant, this understanding of the away-side jet in terms of conical shock waves still 
needs confirmation. A very recent and improved analysis by the STAR collaboration has 
given further support to this picture [l5^. A more solid evidence of this phenomenon may 
come from the study of jets at the Large Hadron Colider (LHC), where the energy released 
by the nuclear projectiles in the central rapidity region will be larger [l6| and so the formed 
fireball will be larger and live longer, allowing for a more complete study of waves. 

In this work we discuss another possible mechanism for the formation of broad structures 
in the away-side jet. In the limit where the jet looses most of its energy, which is rapidly 
thermalized and incorporated to fluid, a pulse is formed, which propagates through the fluid. 
During its motion this energy density pulse spreads both in the longitudinal and transverse 
direction. After hadronization this travelling and expanding "hot spot" will form particles 
with a broader angular distribution than those coming from the near-side jet. This is de- 
picted in Fig. 1. Notice that in this process there is no Mach cone formation. During the 
motion of the energy density pulse, the medium undergoes an expansion leading to a spread 
of this pulse. A further spreading will occur during the hadronization and final particle 
formation. Therefore, in this picture it is essential that the initial perturbation remains 
localized to a good extent. Otherwise it will spread too much and destroy the jet-hke topol- 
ogy, which is compatible with data. Highly localized perturbations can exist and propagate 
through a fluid. The most famous are the Korteweg - de Vries (KdV) solitons, which are 
solutions of the KdV equation. This equation may be derived from the equations of hydro- 
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FIG. 1: Parton-parton collision forming two back-to-back jets, which evolve in a hot quark gluon 
plasma. The circles represent localized (soliton-like) energy density perturbations which traverse 
the fluid and suffer expansion, forming a narrow near-side jet and a broad away-side jet. 



dynamics under certain conditions. One of them is to preserve the non-linear terms of the 
Euler and continuity equations. The other one is to have a third order spatial derivative 
term. This term comes from the equation of state of the fluid and it appears because the 



13, 



m 



19| or because of the Lapla- 



Lagrangian density contains higher derivative couplings 
cians appearing in the equations of motion of the fields of the theory 20!] • This happens, for 
example, in the non-linear Walecka model of nuclear matter at zero and finite temperature. 
For a quark gluon plasma (QGP) it depends on the coupling regime and on the properties of 
the QCD vacuum. As it will be seen in this work, if we consider the simplest case of a free 
gas of massless quarks and gluons, the hydrodynamical equations do not give origin to the 
KdV equation. Instead they generate a non-linear differential equation for the perturbation 
which has no third order stabilising term. This equation is called wave breaking equation 
and is also very well knonw in the literature. The numerical solution of this equation shows 
that an initial gaussian-like perturbation in the energy density evolves creating a vertical 



"wall" in its front, which breaks and looses localization. In our case, surprisingly enough, 
this same phenomenon happens but it takes a very long time and long distances, compared 
to the nuclear scales. So, from the practical point of view, there is no distinction between 
a breaking pulse and a soliton. This persistence of localization in the breaking wave is the 
main result of our paper and gives support to the process shown in Fig. 1. However, from 
this finding to a realistic calculation and a serious attempt to describe the data there is still 
a long way. The next step now will be to quantify the broadening of the moving bubble 
in Fig. 1, which will be directly reflected in the angular distribution of the fragments. For 
this we need to extend our formalism to two spatial dimensions (longitudinal x and radial 
r). This is a heavily numerical project and it is still in progress. Based on previous works 
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HI], very often the hydro- 



with the analogous non-relativistic problem for nuclear matter, discussed in [2l|], we have 
reasons to expect a soliton-like evolution along the x direction with a "leakage" to the radial 
direction, which would cause the angular broadening in the final matter distribution. 

In the theoretical description of these perturbations 
dynamic equations are linearized for simplicity. As it is usually done in non-relativistic 
hydrodynamics, linearization consists |22l] in considering only first order terms in the ve- 
locity and in the energy and pressure perturbations and neglecting higher order terms and 
derivatives involving them. In this work we revisit the relativistic hydrodynamic equations 
expanding them in a different way, in terms of a small expansion parameter (cr) closely fol- 
lowing what is done in magnetohydrodynamics of plasmas 23|] and keeping the non-linear 
features of the problem. Techniques of plasma physics started to be applied to nuclear hy- 
drodynamics long ago [2J, l25| to study perturbations in the cold nucleus, treated as a fluid 



We extended those pioneering studies to relativistic and warm nuclear matter 
and now to the quark gluon plasma (QGP 
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3,Q, 



The most interesting aspect of 17|, 
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24 



25l | was to find at some point of the 



developement, the (KdV) equation for the perturbation in the nuclear matter density. This 
is the "nuclear soliton". Our main contribution was to establish a connection between the 
KdV equation (and the properties of its solitonic solutions) and a modern underlying nuclear 
matter theory (which in our case was a variant of the non-linear Walecka model) and then 
to show that the soliton solution exists even in relativistic hydrodynamics [itI. \m\. 

In the next section we review the main formulas of relativistic hydrodynamics. In section 
III we discuss the quark gluon plasma equation of state. In sections IV and V we show how 
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to derive the diffferential equations which govern the time evolution of perturbations at zero 
and finite temperature respectively. In section VI we present the numerical solutions of the 
obtained differential equations and in section VII we present some conclusions. 



II. RELATIVISTIC FLUID DYNAMICS 

In this section we review the main expressions of one dimensional relativistic hydrody- 
namics. Throughout this work we employ natural units c = 1, h = 1 and (Boltzmann's 
constant) ks = I- The velocity four vector u'^ is defined as m° = 7, m = 'yv, where 7 is the 
Lorentz factor given by 7 = (1 — f^)^^/^ and thus u'^u,y = 1. The velocity field of the matter 
is V = v{t, X, y, z). The energy-momentum tensor is, as usual, given by: 



(1) 



where e and p are the energy density and pressure respectively. Energy-momentum conser- 
vation is ensured by: 

a.T/ = (2) 
The projection of (jj) onto a direction perpendicular to gives the relativistic version of 



the Euler equation 



26 



,^ 1 
dt ^ ^ {e + p)l 



dp 
'dt 



The relativistic version of the continuity equation for the baryon density is 



9.JB^ = 



26|: 



(3) 



(4) 



Since jb'^ = u'^Pb the above equation can be rewritten as: 



(5) 



The relativistic version of the continuity equation for the entropy density is given by the 
projection of ([2]) onto the direction of m'^ [27|: 



{e + p)d^u^' + u^d^e = 



(6) 



At this point we recall the Gibbs relation: 



e +p = hbPb + Ts 



(7) 



and the first law of tliermodynamics: 



de = Tds + HBdpB 



(8) 



We will later consider a hot gas of quarks and gluons, where the net baryon density is zero, 
i.e., p_B = {dpB = 0) at T 7^ 0. Using this last relation in ([8]) and then inserting ([8]) and 
([7]) in ([6]) we arrive at 



which was expected for a perfect fluid. For future use, the above formula will be expanded 
as: 



which is quite similar to ([5]). 

III. THE QGP EQUATION OF STATE 

We shall use a simple equation of state derived from the MIT Bag Model. It describes 
an ideal gas of quarks and gluons and takes into account the effects of confinement through 
the bag constant B. This constant is interpreted as the energy needed to create a bubble 
or bag in the vacuum (in which the noninteracting quarks and gluons are confined) and it 
can be extracted from hadron spectroscopy or from lattice QCD calculations. There is a 
relationship between B and the critical temperature of the quark-hadron transition Tc which 
is determined by assuming that, during the phase transition, the pressure vanishes. 

The baryon density is given by: 



Ts{dX) + Tu^{d^s) = 



and finally at 



9^(su^) = 



(9) 




(10) 




(11) 



where 



1 



(12) 



1 + Q{k~lf^)/T 



and 



^k = ^kiT) 



1 



(13) 



1 + g{k+ltM)/T 



6 



where from now on /i is the baryon chemical potentiaL At zero temperature the expression 
for the baryon density reduces to: 



PB 



3tt' 



(14) 



where fci? is the highest occupied level. The energy density and the pressure are given by: 



and 



e = B 



p=-B 



1 f 7g 



(27 



3 [ (27r) 



/ £k k {e^l^ - ly^ + / S-k k 



(15) 



(16) 



The statistical factors are 7^ = 2 (polarizations) x 8 (colors) = 16 for gluons and 7q = 
2 (spins) X 2 (flavors) x 3 (colors) = 12 for quarks. From the above expressions we derive the 
useful formulas: 



3{p + B) =e-B= — T^ + — dke[n^ + n 



and 



The speed of sound, cs, is given by 



1 4 

p = -e B 

^33 



2 dp 1 
=d-e = 3 



(17) 



(18) 



(19) 



IV. WAVE EQUATION AT ZERO TEMPERATURE 

In the core of a dense star the temperature is close to zero and the baryon density is very 
high. The quark distribution function becomes the step function. Using (fill) in ([T5l) and 
(IISD we find: 

.(p^)=g)'^V/W/^ + S (20) 

and 

MpB) = ^g)'Vw/^-S (21) 

From (fTSll we have Vp = |V£ and also ^ = Combining these expressions with (l20l) 
and (I2T]) we find: 

Vp = ^g)'VVB^/^VpB (22) 
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and 



dp 

Finally, substituting ([20]), (El]), ([22]) and ([23]) into ([3 

iv' - 1) 



PB 



- + (..V). 



dt 

we obtain: 



VpB + V- 



dt 



(23) 



(24) 



which is the relativistic version of Euler equation for the QGP at T =^ . 

Following the same formalism already used for nuclear matter in 
now expand both ([S]) and ^IM in powers of a small parameter a and combine these two 



Q,[i9l,Q 



we will 



equations to find one single differential equation which governs the space-time evolution of 
the perturbation in the baryon density. We write ([5]) and (12^ in one cartesian dimension 
ix] in terms of the dimensionless variables: 



P 



(25) 



— V - — 

Po ' Cs 

where po is an equilibrium (or reference) density, upon which perturbations may be gener- 
ated. Next, we introduce the ^ and r "stretched" coordinates 
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25j: 



cr 



1/2 



[x - est) 
R 



a 



3/2 



est 
R 



After this change of variables we expand ( l25l) as: 



p = 1 + api + a P2 + 



(26) 

(27) 
(28) 



V = avi + cr f 2 + . . . 

Neglecting terms proportional to cr" for n > 3 and organizing the equations as series in 
powers of a, ([5]) and (1241) aquire the form: 



dpi 



dvi 



cr 



dv2 



dp2 dpi 
dr 



dvi dpi 



2 dvi 



and 



J 1 (9pi dvi 



9e 







Pi dpi 



1 gp2 dvi 

3cs2 di 9r ' 9e 9e 3 9e 30^2 9^ 







respectively. In these equations each bracket must vanish independently, i.e. {. . .} = 0. 
From the terms proportional to cr we obtain 05^ = 1/3 and pi = fi, which are then inserted 
into the terms proportional to giving after some algebra: 



dpi 2 dpi 







(29) 



Returning to the x — t space the above equation reads: 

dpi . dpi 2 ^ dpi 

where we have used the notation pi = api, which is a small perturbation in the baryon 
density. The equation (!30l) is the so called breaking wave equation for pi at zero temperature 
in the QGP. 

V. WAVE EQUATION AT FINITE TEMPERATURE 

In the central rapidity region of a typical heavy ion collision at RHIC we have a vanishing 
net baryon number, i.e., pb = 0. The energy is mostly stored in the gluon field, which forms 
the hot and dense medium. We will now aply hydrodynamics to study this medium and 
focus on perturbations in the energy density and their propagation. Following the formalism 
developed in the previous section we will expand and combine the Euler equation given by 
([3]) and the continuity equation for the entropy density given by (ITOl) . 

As Pb = 0, the baryon chemical potential is zero (/i = 0) and so the distribution functions 
given by f|T2|) and (fl3|) are the same, i.e. : = = 1 / (1 + e^^'^). In this case the integral 
in ( ITTI) can be easily performed and we obtain: 

?,{p + B)=s-B=^^^^T' (31) 



Solving the first identity for the pressure and recalling 



28| that s = {dp/dT)v we arrive at: 



s = ^(-B+ -tt'tA = A-n'T' (32) 
9TV 90 / 90 ^ ^ 

The "bag constant" parameter, B, can be replaced by the critical temperature of the 
quark-hadron transition T^. When p = at the phase transition, (|3T1) reduces to: 

B = ^TT^IT.)^ (33) 

Inserting the above equation into the second identity of (!3T|) we have the following expression 
for e(T): 

30 V 3 ; 

From this formula we can define the reference energy density Eq, which is related to a 
reference temperature, Tq, through: 



Solving the second identity of (13T|) for the temperature we obtain: 



T 



30 
37^ 



{e-B) 



1/4 



(35) 



which, inserted into ( l32l) yields: 



30 
37^ 



3/4 



(36) 



Substituing then fl36|) in ffTOl) in the one dimensional case and using fl3Tl) to write {e — B) in 
terms of the temperature we have finally: 



(1-.^) 



90 



de dv 



90v \de 



UStt^T^ J dt dx \ UStt^T^ J dx 
Also from fl^ we have 



/ dv dv^ 
\ dt dx 







e +p 



148 

"go" 



2rr4 



vr 1 



Inserting the above equation into ^ and using Vp = |Ve and also ^ = |f| we find: 

' dv dv^ 



148 



de , de\ 
dt) 



We now rewrite fl371) and fl39|) in dimensionless variables: 



(37) 



(38) 



(39) 



e 



cs 



(40) 



where Eq is the reference energy density given by (134|) . Expanding (140|) in powers of cr we 
have: 



and 



£ = 1 + aei + a £2 + • • • 



V = avi + cr f 2 + . . . 



(41) 



(42) 



Neglecting higher order terms in a and changing variables to the (^ — r) space the equations 
dSZD and ^ become: 



a< 



90 en dei dvi 1 



1487r2T4 9^ 9^ 



90 
1487r2T4 



Sea , Se^i , dei\ dv2 2 ^^i 1 n 



(43) 



and 



0"< 



UStt^T'^cs dvi eo dei 

30 ac" ^ ac" 



>+ 
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a 



30 



dv2 dvi dvi 



eode2 dei 
cs 







(44) 



As before, in tlie above equations eacli bracket must vanisli independently. From the first 
bracket of fH5]) we have: 

90eo 



1487r2r4 

which, inserted into the terms proportional to cr^, yields: 



+ 



90er 



2 dei 







Or ' \U8TT^Ty3 " 
Coming back to the x — t space the above equation becomes: 

90^0 \2 ^ dii 



(45) 



(46) 



Oil Oil ^ . 

'dt + '''' a7 + I 148^ J S'^'^^'d^ ^ 



(47) 



where ii = aei is a small perturbation in the energy density. Equation (H7|) is the breaking 
wave equation for ii in a QGP at finite temperature. For our purposes we will rewrite this 
expression in a slightly different form. Using ( !34l) and the relations deduced in the previous 
section, fHTI) becomes finally: 

(48) 



9t dx 



1 Tc 

An, 



C5. dii 
2^' dx 







where Tq > T^. 



VI. NUMERICAL ANALYSIS AND DISCUSSION 



The equations fl30|) and psj) have the form 



dt dx dx 



which is a particular case of the equation: 



dt dx dx dx^ 



(49) 



(50) 



when B = 0. The last equation is the famous Korteweg - de Vries (KdV) equation, which 

I 

has an analytical soliton solution given by [29]: 



f{x,t) = ^-^^^sech' 



a 



(m - Cs) 



AB 



X — ut) 



(51) 
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where u is an arbitrary supersonic velocity. 

A soliton is a localized pulse which propagates without change in shape. On the other 
hand, the solutions of fl3Up and (HSj) will break, i.e., they will aquire an oscillating behavior 
and will be spread out, loosing localization. Whether or not a given physical system will 
support soliton propagation depends ultimately on its equation of state (in our case, on the 
function e = e^ps) oi e = e{p)). If the EOS takes into account the inhomogeneities in the 
system, the energy density will, in general, be a function of gradients and/or Laplacians. 
When used as input in hydrodynamical equations, these higher order derivatives will, after 
some algebra, lead to the KdV equation. In a hadronic phase, where the degrees of freedom 
are baryons and mesons, we have shown Q, Q, Q, Q that the hydrodynamical equations 
will indeed give origin to the KdV equation. In the present case, for this simple model of 
the quark gluon plasma this was not the case and we could only obtain the breaking wave 
equation. 



A. Zero temperature 

Although the main focus of this work are the perturbations in a hot QGP formed in 
heavy ion collisions, for completeness, we discuss in this subsection the zero temperature 
case, which might be relevant for astrophysics. 

We will present numerical solutions of (1301) with the following initial condition, inspired 
by dSH) 

(52) 

where A and B represent the amplitude and width (of the initial baryon density pulse) 
respectively. In Fig. 2 we show the numerical solution of ( !30l) for A = 0.075 and B = 1 
fm for different times. We can observe the evolution of the initial gaussian-like pulse, the 
formation of a "wall" on the right side. Fig. 3 shows the numerical solution of (|30l) for 
A = 0.35 and B = 1 fm. The time evolution of the pulse is similar to the one found in Fig. 
2 but the "wall" formation and dispersion occurs much earlier. In Fig. 4 we present another 
solution of (130!) for A = 0.075 and B = 0.5 fm. We can see that the initial pulse starts to 
develop small secondary peaks, which are called "radiation" in the literature. Further time 
evolution would increase the strenght of these peaks until the complete loss of localization. 
From these figures we learn how the solution depends on the initial amplitude and width: 



Pi{x, to) = A sech^ 



X 

B 
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FIG. 2: Time evolution of a baryon density pulse at zero temperature. 



it lives longer as a compact pulse for smaller amplitudes and larger widths. Changes in one 
quantity may compensate the changes in the other, creating a very stable moving object. 
In fact, the most striking conclusion to be drawn here is that for a wide variety of choices 
in the initial conditions the solution remains stable and localized for distances much larger 
than the nuclear size. 



B. Finite temperature 

We now turn to the study of the solutions of (l48ll for initial conditions given by (!52l) 
(replacing pi by fi). Now, beside the amplitude and width, the solution will depend also on 
the temperature. When Tq = Tc ( HHi) reduces to: 



dii dsi 2 . dii . . 

which, changing the function from ii to pi is equal to (l30ll . The previous conclusions are 
then extended to the present case. When Tq » Tc (l48l) reduces to: 

dii dii cs . dii 
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FIG. 3: The same as Fig. 2 for a larger amplitude. 



Observing these two formulas we can see that, since Cg = 1/3 is fixed, the only change in the 
differential equation with temperature happens in the numerical coefficient of the last term 
which goes from 0.5 to 0.66. Therefore our results depend very weakly on the temperature. 
A stronger dependence on T would appear if Cs was allowed to change with temperature. 
This would correspond to having a different and more complicated equation of state for the 
quark gluon plasma. 

In Fig. 5 we show the solution of (HHI) with the initial condition given by ( l52l) with 
A = 0.01, 5 = 1 fm and T = 300 MeV. Fig. 6 shows the same as Fig. 5 but with A = 0.1 
and B = 1 fm. As in the zero temperature case, we observe that increasing the initial 
amplitude the breaking process develops earlier. In Fig. 7 we show the same as Fig. 5 but 
with A = 0.01 and B = 0.2 fm. Figs. 8 and 9 show the time evolution of a pulse with the 
same initial amplitude {A = 0.5) and width {B = 1 fm) but different temperatures. Even 
though one temperature is T = 150 MeV (Fig. 8) and the other is T = 300 MeV (Fig. 9) 
we can hardly notice any difference. 
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FIG. 4: The same as Fig. 2 for a smaller width. 



VII. CONCLUSIONS 



We have proposed an alternative explanation for the observed broadening of the away- 
side peak. It is based on the hydrodynamical treatment of energy perturbations. In contrast 
to other approaches we went beyond linearization of the fundamental equations and did not 
neglect the non-linear terms. We used a simple equation of state for the QGP and expanded 
the hydrodynamic equations around equilibrium configurations. The resulting differential 
equations describe the propagation of perturbations in the energy density. We solved them 
numerically and found that localized perturbations can propagate for long distances in the 
plasma. Under certain conditions our solutions mimick the propagation of Korteweg - de 
Vries solitons. However, as said before, from this finding to a realistic calculation and a 
serious attempt to describe the data there is still a long way. The main result found in this 
work, namely, the persistence of soliton-like configurations, is very promising and encourages 
us to extend our formalism to two spatial dimensions. This project is in progress. 



15 









0,010- 


is I I ■ 












■■ 1 : 


i \ 

ii li 




T=300 MeV 


0,008 - 




'1 : 
I' :: i! : 


:: i! 




t=Ofm 






' i! i 






t=5fm 






' ■■ ;i : 


li 1 




t=15fm 


0,006 - 




li ; 

i! i 


■ • 

i! 




1=50 fm 






'\ ■■: !■ i 


11 i! 




t=75 fm 








M ii 




t=100fm 


0,004 - 






il 

i : i i 




t=120fm 












t=150fm 


0,002 - 




, : : i ! ; 

i 1 i 
' ■■ i ! : 


i i 1 
i ■ ' ■ 

1 : : i ! 

' i ' 




0,000 - 


1 




''1 .M. n 1 .-^ 


I 


1 ' 1 ' 1 ' 1 



x(fm) 



FIG. 5: Time evolution of an energy density pulse at T = 300 MeV. 
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